International Journal of Engineering Inventions 

e-ISSN: 2278-7461, p-ISSN: 2319-6491 
Volume 3, Issue 5 (December 2013) PP: 21-23 



Modeling the dynamics of molecular concentration during the 

diffusion procedure 

LingZhen Cao 1 , JiHong Yu 2 

College of life science, Jiangxi Normal University, Nanchang, P.R. China 



Abstract: In this paper, a numerical method is proposed to solving a reaction-diffusion system. The systems 
describe the solid state reaction mechanism of all the reaction components at the preparatory stage of yttrium 
aluminum garnet (YAG) in one dimension. Kansa 's method and forward Euler method are employed to discrete 
the spatial variable and time respectively. 
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I. INTRODUCTION 

In this paper, we employ the Kansa's method to solve the nonlinear partial differential equation system 
defined on the space x £ I c R 1 as following. 

dci _ d 2 ci „ 

^ = Di^-3kc iC2 (1) 
^ = D 2 ^-5kc iC2 (2) 
-^ = D 3 -^+2kc iC2 (3) 

with initial conditions q (x, 0) = q° (x) , x £ I = (a, b) and with zero Neumann boundary conditions on the 

dc- 1 

boundary x = a and x = b, that is, — = 0, i = l,2,3,when t > 0. 

ax lx=a Ux=b 

The mathematical model above is a reaction-diffusion system, which describes the solid state reaction 
mechanism of all the reaction components at the preparatory stage of yttrium aluminium garnet (YAG ). Where 
q = q(x,t) (i=l,2,3) are the concentrations of the i-th reactant respectively, i.e., A1 2 0 3 ,Y 2 0 3 and Y 3 A1 5 0 12 in 
the reaction 

3A1 2 0 3 +5Y 2 0 3 ^2Y 3 A1 5 0 12 , 
of the synthesis at a point x £ (a, b) at time t. D are the diffusion coefficients, that is, we 

consider all the diffusion coefficients coincide. This assumption is general when this sizes of molecules are 
nearly the same. For details, we refer the references [1-7]. 

It is well known that the system (l)-(3) in common three-dimensional case is complicated and so it is 
difficult to implement the numerical modeling. In [2], the authors carried out the solution of the system (l)-(3) 
numerically by using finite difference techniques in one-dimensional case and in two-dimensional case. 
Sysmmetric implicit scheme and alternating direction scheme are employed in the former case and in the latter 
case respectively. Both schemes were solved using stream sweeping method [8]. However, to guarantee the 
convergence of the method of finite difference, it is necessary that the Courant-Friedrichs-Lewy (CFL) 
condition, which says that the 

time step should be less than a certain time in many explicit time-marching computer simulations, must be 
satisfied[9]. 
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In this paper, a meshless method is employed to compute the concentrations of the reactant of A1 2 0 3 ,Y 2 0 3 and 
Y 3 A1 5 0 12 in the reaction 

3A1 2 0 3 +5Y 2 0 3 ^2Y 3 A1 5 0 12 . 

Generally, numerical methods, such as the finite element method (FEM), finite differences (FD) and 
finite volumes (FV), transform a partial differential equation over continuum into a finite set of algebra 
equations. In order to accomplish this tast, it is needed to discretize the domain into a grid, mesh or set of points 
with a fixed connection among them. However, the generation of a well-behaved mesh is not a trivial question 
except in very regular geometries. To compensate the difficulties in traditional mesh-based numerical methods, 
meshless methods are emerged. Moreover, meshless methods are very fit for solving the problems in three 
dimension space or when the domain is complicated, because there are no fixed connectivities among the nodes. 
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Meshless methods have gained much attention in recent years, not only in the mathematics but also in the 
engineering community. Much of the work concerning meshless approximation methods is the interdisciplinary 
between mathematics and engineering. Many traditional numerical methods can either not handle such problems 
at all, or are limited to very special situations. Meshless method is often better suited to cope with the changes in 
the geometry of the domain of interest than classical discretization techniques such as finite difference, finite 
elements or finite volumes. Another obvious advantage of meshless discretizations is their independence from a 
mesh. Mesh generation is still the most time consuming part of any mesh-based numerical simulation. Since 
meshless discretization techniques are based only on a set of independent points, these costs of mesh generation 
are eliminated. Meshless approximation methods can be seen to provide a new generation of numerical tools. 
There are many kind of meshless method, here we make a brief introduction of Kansa's method, for more 
details, we refer to the references [10]. 

A function <D: R s -> R is called radial provided there exists an univariate function 4>: [0, +00) -» R such that 

O(x) = q>(r), where r = ||x||. 
|| . 1 1 usually is the Euclidean norm on R s . It is obvious that for a radial function <1>, <l>(x 1 ) = 0(x 2 ) holds for all 
ll x ill = ll x 2ll ( x i' x 2 G R s )- in standard Gaussian radial basis function (RBFs) approximant, the approximate 
form of function C is written as 
c(x)*Z{L 1 a,<|>(||x-x I ||) > xEV. (4) 

where {a;(i = 1,2, ... , N)} £ R are the unknown coefficients to be determined. 
EfLiOiCKIIx-Xjll) = c(x k ),k = 1,2 N. (5) 

It is to be noted that there is a large class of RBF'J 1 available and one have to choose appropriate one according 
to his need. The most studied RBF including the multiquadric (MQ) cf>(r) = Vr 2 + c 2 , inverse MQ cf) (r) = 

(r 2 + c 2 )~ s , thin plate splines cb(r) = r 2 logi@r)and Gaussian expiipi- (0 ]. For more details about RBF, we 
refer the famous review of 29 interpolation method over scattered data [11]. 

In early 1990s, Kansa made the first attempt to extent RBF interpolation with multiquadrics to the solution of 
PDEs in the strong form collocation formulation [12, 13]. The method, named the Kansa's method, is a domain 
type numerical technique in the sense that the problem is discretized not only on the boundary to satisfy 
boundary conditions but also inside domain to satisfy the governing equation. 

II. NUMERICAL METHOD FOR SOLVING THE PROBLEM 

To solve the systems (l)-(3) with the given initial and boundary conditions when I c R x ,we make 
assumption that all particles are of the same shape of cube shape and its volume is small enough. These are the 
same as in the papers [2, 3] and readers can find more details in them. 
In light of Kansa's method, we write the concentrations q, (1 = 1,2,3) as 

c 1 (x,t)=2]f =1 a!(t)4)(x-x i ). (6) 
Here, the set {Xj,i = 1,2,3, ...,N} is consisted by Halton points [14] in I and is chosen as the central set. x £ 
I U(a, b) are known as the test points. 

For simplicity, we introduce some other notations 

Cf); = CKUX-Xill) , 

<K = <KII X - X illlx=x k , 

^(1) _ flcKllx-Xjll) 

- ax? 
o}' m = aj(x,t m ), 

co T ={t m : t m = mx, m = 0,1, ... , M},Mt = T. 

To solve the non-linear system (3), we propose an alternative iteration algorithm. The scheme is given as 

c.^+lWx.t") = D £N =i a .,m + l qOQ + DikZ N =i a l,m ^ £N =i a 2, m ^ ~ 

or equivalently, 

EJU Jr + % = t[d ziU «!' m ^ + D^lt rt m ti)G3U <*?' m foH+ZlLi O, 

Where D t equals to— 3,— 5,2 whenl = 1,2,3, respectively. In the form of system, it can be written as follows, 



ax 

(i) _ a ct>cilx-x,||) 



ax 

(2) _ fl 2 (j>CI|x-Xil|) 



x=x k 



ax 2 

(2) _ a 2 <j>cilx-xiii) 



'x=x k 
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Zl,m+1 , 

i=l 
N 

Z2,m + 1 . 



i=i 

N 

[D^> m 4\ ( 



3k ( 2^ a i ' m 4>i I af' m + a- ' m ct>i 

\i=l / i=l i=l 

/ N \ N N 

5k [Y £ af' m 4)0] + ^ of'-ch 



i=l 
N 



i=l 
N 



vi=l 
N 



i=l 
N 



i=l 
N 



^ a^ m+1 ^ = T [D ^ a; m 4> 1 (2) + 2k £ o^ m 4>i £ af' m 4)0] + £ a?' m 4>i 



(7) 

When considering the boundary conditions, we denote 

(1) a^Clla-Xjll) 

(a) = Fx 

Then we can obtain the following linear systems. 
A'X 1 = F 1 , (8) 

where A 1 £ r(n+2)x(n+2) with ife entry ^ = ^ when k = 12j ...^^ = ^(a) and A' ki = 4> 1 (1) (b) when 
k = N + 1 and k = N + 2 respectively. 

F> = x[D SJLi «l' m (Z?=i ct>f k 2) ) + D.KCZJL! 4> ik )Q]?li af' m 4> ik )] + Z? =1 
when k = 1,2, ... , N and F 1 = 0 k = N + 1 and k = N + 2. 
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